Electron-lattice relaxation, and soliton structures and their interactions in polyenes. 
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Density matrix renormalisation group calculations of a suitably parametrised model of long 
polyenes (polyacetylene oligomers), which incorporates both long range Coulomb interactions and 
adiabatic lattice relaxation, are presented. The l 3 Bi and 2 1 Af states are found to have a 2-soliton 
and 4-soliton form, respectively, both with large relaxation energies. The 1 X B~ state forms an 
exciton-polaron and has a very small relaxation energy. The relaxed energy of the 2 1 state lies 
below that of the 1 .8,7 state. The soliton/anti-soliton pairs are bound. 

PACS numbers:71.10.F, 71.20.R, 71.35 
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Electronic interactions in polyenes and polyacetylene 
(PA) induce strong spin-density-wave correlations in the 
ground state, resulting in low energy spin-flip (or cova- 
lent) triplet ( 3 B+) excitations. These combine to form 
even- parity (dipole- forbidden) singlet ( 1 ^4g ) excitations. 
Optical (dipole-allowed) transitions to the odd-parity sin- 
glet state CB~) are essentially ionic in character, re- 
sulting in charge transfer from one site to another. In 
the non-interacting limit the 1 3 B+ and 1 X B~ states are 
degenerate, and the 2 1 A+ state always lies higher in en- 
ergy. However, electron correlations can lead to a reversal 
of the energetic ordering of the 1 X B~ and 2*^4+ states. 
Electron-electron correlations in 7r conjugated systems, 
such as PA, are conveniently modelled by the one-band 
Pariser-Parr-Pople (P-P-P) model, which includes long 
range Coulomb interactions. 

Elcctron-phonon interactions in the non-interacting 
limit are described by the SSH model. In the adiabatic 
limit it predicts a wealth of non-linear excitations, in- 
cluding charged/spinless (S ) and neutral/spin 1/2 (S a ) 
solitons. It is the inter-play of both electron-electron and 
electron-phonon interactions in PA which leads to an ex- 
tremely rich variety of excitations. To describe these 
excitations we employ the density matrix renormalisa- 
tion group (DMRG) method to solve the P-P-P-SSH 
model, and utilise the Hellmann-Feynman (H-F) theorem 
to calculate the low-lying excited states and the lattice 
relaxation associated with them. 

Earlier work on the solitonic structure of the low-lying 
excitations include, a renormalisation group calculation 
of the Hubbard-SSH model of up to 16 sites M; a mean- 
field study of the Heisenberg-Peierls model |3j; an exact 
diagonalisation of a 12 site extended Hubbard-SSH model 
|Q| ; and a strong coupling and perturbation calculation of 
the Hubbard-SSH model |). The DMRG method has re- 
cently been used by and Yaron et al. || and Fano et al. |?|] 
to solve the P-P-P model for linear and cyclic polyenes, 
respectively. Jeckelmann || studied the metal-insulator 
transition in doped PA by solving the extended Hubbard- 
SSH with the DMRG method. Likewise, Kuwabara et al. 



II used the DMRG method to study the relative stability 
of bipolarons using the same model. 

The P-P-P-SHH Hamiltonian is defined as 
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where, U = (t + 4>) and f l = | E CT (4+1^ + h - c -) 
is the bond order operator of the ?th. bond. We use 
the Ohno function for the Coulomb interaction: = 



[7/yi + /3r 2 -, where [3 = ([//14.397) 2 and bond lengths 

are in A. The single and double bond lengths used in 
the evaluation of Vij are 1.46 A and 1.35 A, respec- 
tively, and the bond angle is 120°. Various semi-empirical 
parametrisations exist for to and U. We adopt the 
values which are optimal for benzene ]l0| , whose C-C 
bond length of 1.40 A is almost the same as the average 
bond length in PA thin films, i. e., to = 2.539 eV and 
U = 10.06 eV. The dimensionless electron-phonon cou- 
pling constant, A, is defined by A = 2a 2 /irKt , where K 
is the elastic spring constant (estimated to be 46 eV A~ 2 
PI), and a relates the actual distortion of the ith. bond 
from equilibrium, j/j, to A*: yi — Aj/2a. T is chosen 
so that the relaxed ground state of an infinite polymer 
has the same chain length as the unrelaxed state, i. e., 
^^J^ 1 Aj = 0. This ensures that the average hopping 
integral is t n , which is applicable to C-C bond lengths of 
1.40 A. However, the chain length is permitted to change 
for excited states, and for all the states of finite oligomers. 
The remaining parameter, A, is chosen so that the model 
fits the vertical excitation energies of the 1 X B~ and 2 1 A+ 
states of hexatriene in the gas phase jL2|. A choice of 
A = 0.115 and T = 0.602 gives 4.965 eV and 5.212 eV, 
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compared to the experimental values of 4.96 eV and 5.21 
eV, for the l x B~ and 2 1 A+ states, respectively. 

The equilibrium values of the bond length distortion 
are determined by the H-F condition that, 
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H possesses spatial reflection, particle-hole and spin- 
flip symmetries. Symmetrised eigenstates of TL are con- 
structed by an efficient process which makes use of the 
fact that the block symmetry operators commute with 
the density matrix at all stages of the calculation and 
which has been tested by making comparisons with ex- 
axt results HQ. 

The calculation of the relaxed energy of a given state 
for a given chain length is as follows: (1) The eigenstate 
is calculated for an initial choice of {ti} by building up 
the lattice to the target chain size using the infinite lat- 
tice algorithm of the DMRG method. (2) At the target 
chain size the H-F condition (Q) is repeatedly applied 
until the {U} have converged. (3) Using the new values 
of {ti}, steps 1 and 2 are repeated. The procedure is 
succesfully terminated when the energies have converged 
after successive lattice and H-F iterations. It is necessary 
to sweep through the lattice after each set of H-F itera- 
tions to ensure that the electronic states and the lattice 
geometry have converged simultaneously. 

The accuracy of the DMRG implementation has been 
checked in a number of ways. First, the method has 
been compared with exact results in the non-interacting 
(U = 0) limit. The convergence of the ground state en- 
ergy with superblock Hilbert space size (SBHSS) is shown 
in Table | for the N = 102 site system with various lattice 
geometries. The total energy converges to within 0.005 
eV which is sufficient for energy gaps, which are of the 
order of 1 eV, to be resolved to within 1% or better. This 
represents the DMRG at its least accurate, as the addi- 
tion of correlations improves convergence, as can be seen 
in Table O, where we present the DMRG convergence for 
the ground state energy and a number of energy gaps. 

Using the ground state geometry, the vertical energies 
(E v ) of the l 3 B+, l l B~ and 2 1 A+ states are calculated. 
These, as well as the relaxed energies (E°'°), are shown 
in Fig. [I] as functions of 1/N. We first note that the ver- 
tical energies of the 1 X B~ and 2 1 states are very close, 
with a crossing at short chains, and again for long chains. 
In the thermodynamic limit E V (1 1 B~) < E V (2 1 A+). 
This large N crossing has also been observed in the U-V 
dimerised Hubbard model 

The relaxation energy of the lrB~ state is modest (ca. 
0.3 eV) and has not converged (i. e. it is still rapidly 
decreasing) for N = 102. By contrast, the relaxation 
energies of the 1 3 -B+ and 2 1 A+ states are substantial, 
being ca. 0.5 eV and 1.0 eV, respectively, and converge 
rapidly with N. We have also calculated the energy of the 
2 1 A+ state using the relaxed geometry of the l x B~ state. 



This always lies lower than E°~° (l 1 B~), which implies 
that a vertical photo-excitation to the 1 X B~ state will 
decay to the 2 1 A+ state. Finally, the experimental values 
of E°-°CV-B-) and E°-°(2 1 A+) for N = 10 and 14 are 



shown |15|. The 2 1 ^4+ values are in good agreement with 
our calculation. The 1 X B~ values are ca. 0.3 eV lower 
than our predictions. The experimental results for N — 
8-14 have been analysed by Kohler |15|. For the 2 1 A+ 
state the empirical relation E°-°(2 1 A+) = 0.96 + 20.72/iV 
was derived, in good agreement with the photoinduced 
absorption result ca. 1.1 eV for polyacetylene thin films. 
However, Fig. 1 suggests that an algebraic fitting form is 
incorrect — the true scaling behaviour is exponential, and 
can only be seen by considering sufficiently large systems. 

In Fig. H we plot as a function of bond index from 
the center of the chain, the normalised staggered bond 
dimerisation, defined as, Si = (—iy(ti — t)/t, where t is 
the average value of t{ in the middle of the chain [fl6|| . 
Note that the 1 3 B+ and 2 1 ^4+ states undergo consider- 
able bond distortion, whereas the 1 X B~ state and the 
charged state (denoted E g ) show a weak polaronic dis- 
tortion of the lattice. The oscillatory behaviour of Si in 
the polaronic distortions indicates a local expansion of 
the lattice. 

We fit the 1 3 B+ , 1 X B~ and charged state to a 2-soliton 
form |@, 
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The 2 1 A+ state, however, evidently requires a 4-soliton 
UfjU fit of the form, 

Si = S{1 + tanh (2xq/£) [tanh {{i — x d — £o)/0 

- tanh ((i ~ x d + #o)/£) + tanh ((i + x d ~ Xq)/£) 

- tanh((i + ar d + ar )/£)]}. (4) 

These functions give good fits to the relaxed geometries 
of the N — 102 site system, as shown in Fig. ||. The 
difference in energy between using the fits and the actual 
relaxed geometry is around 0.01 eV. The 4-soliton charac- 
ter of 2 1 j4+ state indicates the strong inter-play between 
electron-lattice relaxation and electron-electron correla- 
tions in polyenes, for, as indicated earlier, this state has 
a considerable triplet-triplet character. 

Fig. H depicts the convergence of the various fitting pa- 
rameters as a function of N. The 1 3 -B,t and 2 1 ^4+ states 
converge rapidly with N, whereas the 1 B~ state shows 
strong finite-size effects, and the coherence length £ only 
begins to converge at around TV = 102. The fact that 
the soliton structures converge with TV leads to two im- 
portant observations: First, the soliton/anti-soliton pair 
are bound, because if they were not their separation xq 
would increase with N. Second, the soliton structures 
are pinned in the middle of the lattice. This is a conse- 
quence of the classical adiabatic treatment of the lattice, 
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and is one of the reasons why the energy curves flatten 
off rapidly as N — > oo. 

To further investigate the soliton/anti-soliton interac- 
tions, adiabatic potential energy curves Q (i. e. the en- 
ergy as a function of soliton separation, Xq) are plotted 
in Fig. ^. Our results differ qualitatively from previous 
approximate calculations ]Bj in that the ^B~ and l 3 B+ 
are bound — the potentials have a minimum and are at- 
tractive for large xq. This attractive soliton/anti-soliton 
interaction implies much stronger binding for the 
state than that obtained in @] , where the binding energy 
was found to be ca. 0.05 eV. It should be noted, however, 
that § uses a Hamiltonian with short ranged (on-site) 
interactions, with the strength chosen so as to fit the 
vertical absorption peak in polyacetylene thin films. Fur- 
thermore, our calculations, being for polyenes, necessar- 
ily use open boundary conditions. For an even site chain 
there is one more short bond than there are long bonds. 
This means that the ground state is non-degenerate, as 
it is energetically unfavourable to swap long and short 
bonds, and is one reason for the long range confinement. 
The role played by boundary conditions is subtle and im- 
portant, as real systems, such as oligomers and polymers 
with disorder, have a finite conjugation length. 

Another consideration is the adequacy of the soliton 
fits used in generating the adiabatic energy curves. We 
consider the generalised potential energy curves where, 
for a given xq, we allow £ to vary so as to minimise the 
energy. Results for the ^B~ state are included in Fig. ||, 
which show that relaxing £ yields a substantial reduction 
in the energy, implying weak soliton/anti-soliton bind- 
ing. However, the energy reduction for the 1 3 -B^~ state is 
insignificant over the range of xq values plotted, indicat- 
ing that there is a stronger binding of the solitons in the 
2 1 A+ and 1 3 B+ states. A further generalisation of the 
soliton fits would be to consider multiple soliton/anti- 
soliton pairs. 

Finally, we note the consequences of our results for 
the interpretation of experiments. Our results for small 
polyenes are in good agreement with experiment — the en- 
ergy difference of ca. 0.3 eV for the ^Bz state can prob- 
ably be explained by solvation effects [|19[ , supporting the 
notion that the covalent 2 1 A+ state is less polarised than 
the ionic l 1 ^-. In the bulk limit the l l B~ and 2 1 A+ 
energies are ca. 0.8 eV higher than data from linear |2C| ] 
and 2-photon pl|| absorption and third harmonic genera- 
tion [p2| experiments on PA thin films. This implies that 
there are more substantial energy decreases due to sol- 
vation and aggregation (interchain hopping and eximcr 
formation) effects. Such effects must be investigated via 
coupled chain calculations. Also, the neglect of quantum 
fluctuations in the adiabatic treatment of the lattice [^3) , 
leading to the pinning of the soliton structures, will con- 
tribute to this energy difference. A full treatment must 
include dynamical phonons. Such a treatment would also 
increase our understanding of the soliton confinement. 



R. J. B. was supported by the Australian Research 
Council. W. B. gratefully acknowledges financial support 
from the Royal Society and the Gordon Godfrey Bequest 
of the UNSW. We thank E. Jeckelmann and Y. Shimoi 
for useful discussions. The calculations were performed 
at the New South Wales Center for Parallel Computing. 



* Email address: phlrb@newt.phys.unsw.edu.au; 

w.barford@shefneld.ac.uk 
[1] S. R. White, Phys. Rev. Lett. 69, 2863 (1992); Phys. Rev. 

B 48, 10345 (1993). 
[2] G. W. Hayden and E. J. Mele, Phys. Rev. B 34, 5484 

(1986). 

[3] J. Takimoto and M. Sasai, Phys. Rev. B 39, 8511 (1989). 
[4] J. T. Gammel and D. K. Campbell, Synth. Met. 55, 4638 
(1993). 

[5] W. P. Su, Phys. Rev. Lett. 74, 1167 (1995). 

[6] D. Yaron, E. E. Moore, Z. Shuai, J. J. Bredas, J. Chem. 

Phys. 108, 7451 (1998). 
[7] G. Fano, F. Ortolani, L. Ziosi, J. Chem. Phys. 108, 9246 

(1998). 

[8] J. Jeckelmann, Phys. Rev. B 57, 11 838 (1998). 
[9] M. Kuwabara, Y. Shimoi, S. Abe, J. Phys. Soc. Jpn. 67, 
1521 (1998). 

[10] R. J. Bursill, C. Castleton, W. Barford, Chem. Phys. Lett. 

294, 305 (1998). 
[11] E. Ehrenfreund, Z. Vardeny, O. Barfman, B. Horovitz, 

Phys. Rev. B 36, 1535 (1987). 
[12] W. M. Flicker, O. A. Mosher and A. Kuppermann, Chem. 

Phys. Lett. 45, 492 (1977). 
[13] M. Boman, R. J. Bursill and W. Barford, Synth. Met. 85, 

1059 (1997). 

[14] M. Boman and R. J. Bursill, Phys. Rev. B 57, 15 167 
(1998). 

[15] B. E. Kohler, J. Chem. Phys. 88, 2788 (1988). 

[16] We have checked the DMRG convergence of the Si and 

found that, for the N — 102 site lattice, they are resolved 

to within 1% or less. 
[17] Z. Shuai, J. L. Bredas, S. K. Pati and S. Ramasesha, Phys. 

Rev. B 56, 9298 (1997). 
[18] D. K. Campbell and A. R. Bishop, Nucl. Phys. B 200, 297 

(1982). 

[19] E. Moore, B. Gherman B. and D. Yaron, J. Chem. Phys. 

106, 4216 (1997). 
[20] Relaxation in Polymers (p 174), ed. by T. Kobayashi. 

World Scientific (Singapore) 1993. 
[21] C. Halvorson and A. J. Heeger, Chem. Phys. Lett. 216, 

488 (1993). 

[22] W. S. Fann, et al, Phys. Rev. Lett. 62, 1492 (1989). 
[23] R. H. McKenzie and J. W. Wilkins, Phys. Rev. Lett. 69 
1085 (1992). 



3 



TABLE I. Convergence of the ground state energy (in eV) 
as a function of the SBHSS for the N = 102 site system in the 
non-interacting case (U = 0) for three geometries defined by 
the soliton form (^), taking £ = 4.03 and (i) xq = (uniformly 
dimerised geometry), (ii) xo — oo (uniformly dimerised with 
long and short bonds reversed), and (iii) xo = 24.87 (a ge- 
ometry with a kink/anti-kink pair placed 1/4 and 3/4 of the 
way along the lattice) . m is the number of states retained per 
block. 



m 


SBHSS 


xo = 


Xo = oo 


x = 24.87 


75 


5920 


-332.57146 


-330.205 


-331.333 


100 


9384 


-332.57217 


-330.403 


-331.366 


150 


22392 


-332.57257 


-330.426 


-331.407 


200 


37512 


-332.57271 


-330.439 


-331.422 


230 


52312 


-332.57272 


-330.446 


-331.428 


270 


72392 


-332.57273 


-330.448 


-331.430 


EXACT 




-332.57276 


-330.452 


-331.434 



TABLE II. Convergence of the ground state (l 1 At) energy 
and vertical and 0-0 transition energies of the 2 1 and 1 
states as a function of the SBHSS for the N = 102 site system. 



SBHSS 



lMj 



2M+ 2M+(0-0) l 1 ^- ^-(O-O) 



15844 
25492 
36312 
54916 
67240 



-509.6330807 
-509.6330971 
-509.6331002 
-509.6331009 
-509.6331010 



2.8927 
2.8795 
2.8764 
2.8744 
2.8737 



1.8051 
1.8008 
1.7972 
1.7963 
1.7959 



2.7719 2.6785 

2.7650 2.6483 

2.7617 2.6392 

2.7605 2.6345 

2.7601 2.6336 
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FIG. 1. Energy gaps for the l 1 ^, 



0.1 



0.125 



0.15 



(solid lines) 2 1 A+ (dot- 
ted lines) and l^Bt (dashed lines) states as a function of 1/N. 
Vertical/0-0 transitions are indicated by thin/thick lines. Ex- 
perimental 0-0 energies of the 1^B~ (diamonds) and 2 1 A£ 
(triangles) states for polyenes in hydrocarbon solution [ Jl5[ , 
The empirical fitting form £ M (2M+) = 0.96 + 20.72/iV, de- 
rived in |Q from the 8-14-site oligomer data is also plotted 
(dot-dashed line). 
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FIG. 2. The geometries [5% as a function of bond index i 
from the center of the lattice) of various excitations: l 1 B^ 
(open diamonds), l 3 Bt (filled diamonds), 2 1 A~t (stars) and 
the charged state, E g , (open triangles), for the N = 102 
site system. The solid lines are fits to the 2-soliton form ([]) 
(and the 4-soliton form (^) for the 2 1 A~t). The inset shows 
the two-point averages ((2i + l)/2, (<5, + <$i+i)/2) for the po- 
laronic E g and 1 B~ states, which are well described by the 
2-soliton fits. 
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FIG. 3. The convergence of the soliton fitting parameters 
xo and £ for the 1 X B~ and 1 3 -B„ states with the lattice size, 
N. The inset shows Xq, £ and Xd for the 2 i A^ state. 
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FIG. 4. Potential energy curves (solid lines) for the l 1 B~ 
(diamonds), 2 1 (stars) and l s Bi (triangles) states for the 
N = 102 site lattice. The dashed curves are the correspond- 
ing ground state (l 1 A^) potential energies. For the l 1 73" 
and 1 3 B+ cases, the curves are generated using the soliton 
pair form ^ with the fitted values of £ (4.01 and 12.12, re- 
spectively) and varying xq. For the 2 1 Ag case the curves are 
generated using the 4-soliton form (^) with the fitted value 
£ = 4.93 and varying xo- Xd is chosen so that the ratio Xd/xo 
remains fixed at its fitted value of 1.35. The solid diamonds 
are the values of the ^B~ energy when £ is also allowed to 
vary. The energies of the vertical, 0-0 and emission transi- 
tions, and the relaxation energy can be read off from this 
plot. This is illustrated using arrowed vertical lines for the 
2 1 A+ state. 
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